Bremsstrahlung in dark matter annihilation 



Vernon Barger^, Wai-Yee Keung^, Danny Marfatia^'^ 

^ Department of Physics, University of Wisconsin, Madison, WI 53706, U.S.A. 
^Department of Physics, University of Illinois, Chicago, IL 60607, U.S.A. 
^Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, U.S.A. 



Abstract 

We show that the energy spectra from dark matter (DM) annihilation into /^/"7 (where 
/ is a hght fermion) via chirahty preserving interactions are identical in 3 scenarios: (a) DM 
is a Majorana fermion and the particle exchanged is a scalar, (b) DM is a Majorana fermion 
and the particle exchanged is a vector and (c) DM is a scalar and the particle exchanged is a 
fermion. For cases (a) and (c), we also calculate the differential cross section to £^f~V, where 
V = W, Z, 7, and £ is a light fermion that may or not be the same as /. The form of the cross 
section depends on whether the DM is a Majorana fermion or a scalar but in both cases its form 
is independent of V. 
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Intensive searches for dark matter (DM) via their annihilation signatures [T] have spurred in- 
terest in the effects of electromagnetic and electroweak bremsstrahlung on the rates for particular 
final states and on the resulting spectra [21 H] . 

In this Letter, we show that for 3 different interactions that are chirality preserving and appear 
in popular scenarios, dark matter annihilation to /^/~7 in the static limit gives identical energy 
spectra. We then connect the common structure for the amplitude to an operator basis. Finally, 
for 2 of the 3 cases, we derive the differential cross section to £^ f~V, where V = W,Z,j and find 
its form to be different for Majorana DM and scalar DM, but independent of V. 

Electromagnetic bremsstrahlung 

We consider the annihilation final state e~^e~j for the interaction terms in the table below. The 
exchanged particles S, W and E (collectively denoted by X) are taken to be negatively charged 
and heavier than the corresponding DM particle D = and (j). We assume that D \s a. standard 
model singlet and its coupling to the electron to be chiral with the interaction involving cr. For 
e^, straightforward replacements can be made. 





Model 1 


Model 2 


Model 3 


DM 


Majorana x 


Majorana 


Scalar (j) 


Exchanged particle X 


Scalar 5 


Vector W' 


Fermion E 


Interaction £/ 


g'S^XdR 




9'E(t)eR 



The super symmetric case of neutralino DM with selectron exchange is a good example of 
Model 1. The case of a new heavy charged W' gauge boson that couples cr and a right-handed 
heavy Majorana neutrino N in variations of left-right theory falls into the class of Model 2 {5]. 
Scenarios of scalar DM which have been of interest recently are realizations of Model 3 [1]. 

For the annihilation process, DD — ?• e{pi)e{p2)^{k,e), in the nonrelativistic limit we find a 
unified description of the energy spectra of the decay products for all three models. In standard 
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notation, the amplitude turns out to be 

[{Pi -P2 + - AmW [{pi - P2- ky - 4:m\] ' 

where the minus sign applies to Models 1 and 2 (with Majorana DM) and the plus sign applies 
to Model 3 (with scalar DM). Note that for the Majorana DM models only the initial state with 
vanishing total angular momentum participates in the annihilation. The overall constant 



C 



-^^'2 for Model 1, 

A^g'\2 + ^) for Model 2, 
c/'2 for Model 3. 



Remarkably, we find the same structure for all 3 Models even though Model 2 involves the tri- 
gauge boson vertex W'W'^. The amplitude obeys QED gauge invariance and its two terms do 
not interfere when photon polarizations are summed in the limit mg — )• 0. The spin-averaged 
annihilation rate is 

da ^ 1 |eCp [{I - xrf + {I - X2?]{1 - X,) 

'"^dxidxg (2Sd + 1)2 47r3m2, (l-2xi-r)2(l-2x2-r)2 ' 

where r = m\/m?j^. The scaling variables xi = Ei/rriD [i = 1,2), x^ = E^/m£) are defined in the 
static center of mass frame so that xi + X2 + X3 = 2. The spin-averaged factor (25^+1)^ 3 
Models 1, 2 and unity for Model 3. The photon energy distribution is obtained by integrating over 
XI G [1-X3,1] [2J: 

da 1 |eC|2 I-X3 



dx3 (25d + 1)2 327r3m2, (1 + r - xg)^ 

X f2x3;^ + (l+--^3)^^ _ i^+r)il+r-2x3) 1+r \ 

\ ''(l + r)(l+r-2x3) l + r-X3 1 + r - 2x3 7 

The unified formulas for the amplitude and the distribution for the 3 models are simple and 

interesting. We now study the structure of the amplitude. The numerator of the amplitude involves 

two pieces, each of which are QED gauge invariant. The first one is 

1 _ 

2^«bi) h A /=vr{p2) = ur{pi){p2 -e /:-p2-k ji)vR{p2) < — O , 

where we have used the massless fermion on-shell condition and the last step identifies the stucture 
with an operator [1], 

O = Ff^'^ijR-f.d^lPR . 

Similarly, the second piece is 

^ur{pi) /k / i>iVR{p2) = ur{pi){pi -e /o-pi-k ^)vr{p2) i — , 
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where [Ij 

is the conjugate of O. On the other hand, we can use the Chisholm identity 
to write 

ur{pi) i)2 / /kvR{p2) = ur{pi){p2 ■€/t-p2-k^ + k-€^2 -^e^^'^'^T^Tsj^^fife) 

= -2ieP''''''uRipih^VRip2) . 

Similarly, 

ur{pi) k ^ hvR{P2) = +'2ie^'^^''''uR{pi)jf,VR{p2) . 
So we have the relations, 

O = F^^'^i^R^.d^i^R = -iFf^^i^R-f^df^lPR , 

Ot = F^"'{^f,^PR)-f,^PR = iF'"' id^i;R)^,i;R , 
where F'"' = ie^'^^^Fa/?. Then, 

O + O^ = F^'''d^{ijRjM=T , 

i{0 - Ot) = F'^^df.i'iPRjM = T . 

Thus, either O and O^, or and T can be used as bases to describe the amplitude. Because 
the dual substitution interchanges E o B, it also interchanges a right-polarized photon for a left- 
polarized photon. If the energy distribution sums up polarizations, both F and F give the same 
result. Note that since T is CP even, it corresponds to Model 3 with scalar DM. On the contrary, 
the Majorana pair is CP odd, and the amplitude picks up T . 

Electroweak bremsstrahlung 

For the annihilation of Majorana DM % via scalar exchange, we derive a univeral amplitude for 
vector boson V emission X + X ~^ /(Pi)^(P2)^(fc, e) in the static limit. The vector boson can be 
py, Z or 7. As Vrei 0, only the composite of XX with zero total angular momentum gives a 
contribution, and we get an amplitude of the form of Eq. ([T|) with the minus sign selected, and C 
replaced by Cy (in which we collect the model-dependent couplings and coefficients). As usual, we 
assume that the chirality of the massless lepton field is conserved. Note that the two terms in the 
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numerator do not interfere only when k'^ = (as for the photon V = QED-Uke gauge invariance 
(i.e. • A: = 0) is satisfied even for k"^ = niy 7^ 0. We calculate the annihilation cross section to be 

^ ^ da ^ leCyp [(1 - XI - Syf + (1 - - 5v)^ - 2^(1 - X3 + - X3 + Sy) 

dxidxs IQ-K^m'^ (1 — 2xi — r)2(l — 2x2 — ' 

where r = m^/m^ and 6y = ^m'y/m'^. In the case of W/Z emission, we assume a common mass 
mx for all SU{2) partners of the heavy exchanged particles. Note that 3:3 = Ey/m^ and xi lie in 
the ranges [2y/5y, l + 5y] and [0, 1 — 5y], respectively. The energy distribution of V can be obtained 
by integrating over xi within the limits = (2 — X3 it \J x\ — 45y)/2, and the energy distribution 
of / can be obtained by integrating over X3 from 1 — xi + 5y/(l — xi) to 1 + • On setting 5^=0 
we recover the cross section for V = ^ and is simply C given in the photon case. The values of 
Cw and Cz are given by 

I^W — W n ■ 2 /I ' V'^/ 



ri2 _ ^2 v-^a - \! ^^^^ ^wj fr^-s 



and 

sin^ 9w cos^ 9^/ 

where T, T3, T3 are the weak-isospin numbers of the exchanged particles and 9w is the weak mixing 
angle. Since the spinors in Eq. ([T|) can have either L or it! chirality depending on the underlying 
physics, the weak-isospin numbers must be chosen accordingly. The exchanged particle 5 of Model 1 
is an SU{2) singlet with Q = —1. 

For scalar DM <j) that annihilates via fermion exchange, a relative sign is flipped in the numerator 
of the amplitude, and 2 additional terms appear for V = W,Z, but not for V = 'y: 

^ 8eCyUh{pi){h ^ h+ h ^ h+ ^ /t Jk-k"^ ^)vh{P2) 
[{pi -P2 + kf - 4m^] [{pi -p2- ky - \m\\ 

where h = L,R, and Cw and Cz are as in Eqs. ([5]) and ([6]). The corresponding annihilation cross 
section is 

da ^ \eCy\^ N 

''■■"^dxidx3 47r3m2 (1-2x1 -r)2(l- 2x2- r)2 ' ^ ^ 

iV = [(1 - Xi - 6yf + (1 - X2 - Syf + 26y{l + X3 - 2V)](1 - ^3 + 6y) 

+ 26y(l - Xi - 6y){l - X2 - dy) , 

with kinematic variables defined as before (and replaced by m^). For the interaction Lagrangian 
of Model 3 in the table, the exchanged particle E is an SU (2) singlet with Q = —1. The realization 
of Model 3 in Ref. [4J uses an 5'C/(2)-doublet {N^,E^) for the exchanged particles with T = i, 
^3 = -^, T;^ = ^andQ = -l. 

The difference in the W/Z spectra for Models 1 and 3 arises from terms proportional to 5y. 
In both models, the energy distributions for the photon and the W/Z differ substantially due to 
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kinematic effects if D is just above the W/Z thresholds, while for a high mass D the difference is 
suppressed by 5v- 
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